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Supersymmetric flavor models for the radiative generation of fermion masses offer an alternative 
way to solve the SUSY-CP problem. We assume that the supersymmetric theory is flavor and CP 
conserving. CP violating phases are associated to the vacuum expectation values of flavor violating 
susy-breaking fields. As a consequence, phases appear at tree level only in the soft supersymmetry 
breaking matrices. Using a U (2) flavor model as an example we show that it is possible to generate 
radiatively the first and second generation of quark masses and mixings as well as the CKM CP 
phase. The one-loop supersymmetric contributions to EDMs are automatically zero since all the 
relevant parameters in the lagrangian are flavor conserving and as a consequence real. The size of 
the flavor and CP mixing in the susy breaking sector is mostly determined by the fermion mass 
ratios and CKM elements. We calculate the contributions to e, e' and to the CP asymmetries in the 
B decays to ipK s , <f>K s ,ri' K a and X s j. We analyze a case study with maximal predictivity in the 
fermion sector. For this worst case scenario the measurements of Amjf ,Ams and e constrain the 
model requiring extremely heavy squark spectra. 

PACS numbers: 



I. INTRODUCTION 

It was suggested a few years ago that the huge number 
of possible string theory vacua in combination with eter- 
nal inflation may allow us to understand the smallness 
of the cosmological constant from an anthropic point of 
view [10. Although there is no current general frame- 
work for examining these metastable vacua in string the- 
ory some particular methods have been proposed Q. 
It is still an open debate whether the landscape does or 
does not predict high scale supersymmetry (SUSY) 0,0- 
Statistical analysis of the vacuum of certain string theo- 
ries have derived formulae for the distribution of vacua, 
which support the idea of a very high scale of supersym- 
metry breaking 0. It has been pointed out that this 
considerations may also be relevant to understand the 
gauge hierarchy problem Q . This has motivated the re- 
cent interest in field-theoretic realizations of models with 
large numbers of vacua ||| as well as in the analysis of the 
collider, cosmological and other phenomenological impli- 
cations of supersymmetric models with very heavy su- 
persymmetric spectra [T(ij |. 

In this paper we would like to revisit, under the light 
of these new considerations, certain supersymmetric fla- 
vor models that are not usually considered in the liter- 
ature because it is naively expected that they require 
a very heavy supersymmetric spectra to be compatible 
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with experimental constraints on flavor changing pro- 
cesses (FCNC). In particular we will show that super- 
symmetric flavor models for the radiative generation of 
fermion masses generate very predictive Yukawa textures 
and at the same time offer a new insight in the SUSY CP 
and flavor problems. 

After more than 30 years of its observation, the viola- 
tion of CP symmetry is poorly known in its origin in the 
present particle physics paradigm, despite its relevance 
in nature. The presence of CP violating phases in par- 
ticle physics models can be tested, for instance, through 
precision measurements of the electric dipole moments 
(EDMs) in the leptonic sector (electron and muon) and 
in the quark sector (neutron and deuteron). At present 
there are very stringent upper limits ^lj, which are ex- 
pected to improve by several orders of magnitude in the 
near future experiments. It is known that the Standard 
Model contribution to the neutron EDM, in the absence 
of a theta term, is approximately 10 -30 e cm, which is 
still more than four orders of magnitude below the reach 
of the current experiments. 

In the context of an unconstrained minimal supersym- 
metric standard model (MSSM) 12] the generic contri- 
bution to the electric dipole moments is several orders 
of magnitude larger than the SM contribution. This se- 
rious volation of the current experimental constraints is 
known as the SUSY CP problem. Several possible expla- 
nations have been considered in the literature to account 
for the suppression of the supersymmetric contributions 
to EDMs and other CP violating observables. Some of 
them are: 1) CP suppression 2) cancellations, 3) align- 
ment, 4) sfermion decoupling, and 5) flavor off-diagonal 
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CP violation. 

• The CP suppressed scenario assumes that all the 
CP phases are suppressed p| because CP is an 
approximate symmetry of the full theory. 

• The cancellation scenario is based on the exis- 
tence of certain regions of the SUSY parameter 
space where different contributions to EDMs cancel 

dm 

These two possibilities are nowadays ruled out. The CP 
suppressed scenario would imply that all SM contribu- 
tions to CP asymmetries are small, which we know today 
not to be the case in the B system [T(J . The cancella- 
tion scenario is also known to be ruled out if constraints 
from electron, neutron and mercury atom EDMs are im- 
possed simultaneously 0] . Moreover there seems to be 
no symmetry that can guarantee such cancellations. 

• In the CP alignment case the phases associated 
to the relevant parameters are somehow related in 
such a way that the combinations contributing to 
EDMs cancel. 

• Decoupling entails that the sfermion masses are 
heavy enough to strongly suppress the supersym- 
metric contributions even tough the CP phases can 
be arbitrarily large 01 • 

The alignment scenario could arise naturally in the con- 
text of models that generate exact soft universality as 
gauge mediated supersymmetry breaking models |l9j or 
in models with approximate horizontal abelian flavor 
symmetries |20j . The decoupling scenario is very plausi- 
ble. It requires sfermion masses of the order of several 
TeV, which implies the existence of fine tuning in the soft 
supersymmetry breaking sector. 

• The scenario with flavor off diagonal CP violation 
assumes that the origin of CP violation is closely 
related to the origin of flavor structures in such a 
way that the flavor blind quantities as the ^t-term, 
soft bilinear terms gaugino masses are real and only 
flavor off-diagonal CP phases are non-zero. 

We find that the scenarios with flavor off-diagonal CP 
phases are especially interesting. The models of this kind 
proposed in the literature to date 0, H3| assumed 
that flavor violating Yukawa matrices and soft terms are 
both generated simultaneously at tree level at very high 
energies. Thus, they require that Yukawa matrices and 
soft trilinear matrices are hermitian, which forces flavor- 
diagonal phases to vanish (up to small RGE corrections) . 
We would like to propose a new idea similar to the fla- 
vor off-diagonal scenario which has not been considered 
before, 

• We propose that the underlying supersymmetric 
theory is exactly CP conserving while CP phases 
are only carried by flavor violating susy breaking 
fields. 



At first sight one may be tempted to think that this sce- 
nario cannot account for the observed CP violation in 
the SM, especially the large mixing in the B-B system, 
and therefore conclude that CP violation must be present 
in the superpotential. We will show in this paper that 
this is not the case and certain models for the radiative 
generation of first and second generation fermion masses 
and mixings recently proposed 0, H3 allow us to gener- 
ate radiatively the CKM phase and offer an alternative 
solution for the SUSY CP problem. 

In this model flavor and CP violation appear at 
tree level only in the soft supersymmetry breaking pa- 
rameters and are transmitted to the fermion sector at 
one loop through low energy finite threshold correc- 
tions. CP violating phases could appear originally in 
the vacuum expectation values of certain flavor violat- 
ing susy breaking fields. These vevs break spontaneously 
both flavor and the CP symmetry generating at tree 
level flavor and CP violating soft mass matrices. This 
class of models make use of the presence of soft super- 
symmetry breaking terms for the radiative generation 
of quark and lepton masses through sfermion-gaugino 
loops, as originally suggested by W. Buchmuller and 
D. Wyle r [2a. 1271 a nd later analyzed in more detail in 
Refs. |Sl2lMl3l|3l|3l|3|. The gaugino mass would 
provide the violation of fermionic chirality required by a 
fermion mass while the soft breaking terms provide the 
violation of flavor and CP symmetries. 

In this paper we have chosen to analyze a case study 
that as we will show is the worst case scenario from the 
point of view of FCNC constraints. We analyze a model 
of this kind because it achieves maximal predictivity in 
the quark sector. This case study, as we will see, tends to 
generate important contributions to some flavor changing 
processes, especially Am^, Ato^ and e, which can only 
be avoided if both the squark and the gluino spectra are 
very heavy. 

This paper is organized as follows. We begin in Scc.lTTI 
by describing the model we propose for the radiative gen- 
eration of first and second generation fermion masses and 
mixing angles. In Sec. IIIII we analyze the radiative gen- 
eration of Yukawa couplings in this model. In Sec. IIII Al 
we study the predictions and constraints for quark mass 
ratios. In Sec. IIII Bl we study the radiative generation of 
the SM CKM phase and the predictions and constraints 
arising from measured CKM elements and CP phases. In 
Sec0we argue that the contributions to EDMs are ex- 
actly zero in this model. In Sec. II VI we analyze in detail 
the calculation of the soft matrices in the SuperCKM ba- 
sis. In Sec. I VII we study the contributions to direct and 
indirect CP violation in the K-K system. In Sec. IV III we 
study the contributions to CP asymmetries in the B-B 
system. 
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II. THE MODEL 

In this section we will consider a realistic three gener- 
ation model based in a horizontal XJ(2) H symmetry 13^. 
This is a generalization of the model proposed in Ref. |24| . 
We will assume the usual MSSM particle content where 
third generation matter superfields, 



Q 3 ,V 3 ,U 3 ,£ 3 ,£ 3 



(1) 



and up and down electroweak Higgs superfields, H. u and 
Hd, are singlets under U(2) H . We will denote them ab- 
breviately by 4> L and 4> R - Let us assume that first and 
second generation left handed superfields, 



*o = 



Qi 

Q.2 



C-2 



(2) 



as well as the first and second generation right handed 
superfields, 



Ui 
U 2 



Si 

£2 



(3) 



transform as covariant vectors under \J(2) H , We will de- 
note them abbreviately by ^ and We will introduce 
a set of supersymmetry breaking chiral superfields, 



gab 



A ab , T a (a, b = 1,2), 



(4) 



that transform under XJ(2) H contravariantly as a sym- 
metric tensor, an antisymmetric tensor and a vector re- 
spectively. We will assume that only the auxiliary com- 
ponents of the flavor breaking superfields are non zero, 
The most general form for the vevs of the flavor breaking 
fields is, 



<S> = 
< A > = 
<T> = 









V s e^ s 



U v a e 

-v a e zrt >" 



(5) 
(6) 
(7) 



where the v s , Vs, v a , Vf and Vjr are real parameters. 
We will assume the following particular hierarchy in the 
flavor breaking vevs, 



(v f ,v a ,V s ,Vj:) = (A 2 ,A 2 ,A,2A)M F m. 



(8) 



We will also assume that v s <C Vs and for practical pur- 
poses set v s = 0. Here A is a flavor breaking perturbation 
parameter, M-p is the flavor breaking scale. We note that 
m is a new mass scale linked to the flavor violating susy- 
breaking fields. We do not have yet a predictive model 
for the \J(2) H breaking, which is a relevant point under 
current investigation. The proposed vevs in Eq. [8] are 
introduced ad-hoc. These ad-hoc assumptions will prove 
a posteriori to be very succesfull in reproducing fermion 



masses and mixings. Furthermore, in the case that \J(2) H 
is a gauge symmetry broken spontaneously we expect the 
U(2)^-gauge fields to get masses of the order of the flavor 
breaking scale which can be very heavy in this scenario. 
Therefore any other phenomenological effects in the low 
energy model beyond the flavor structure it gives rise to 
in the soft supersymmetry breaking sector would be very 
suppressed. 

We will assume that the superpotential of the model is 
CP and U(2) H symmetric. Therefore the only couplings 
allowed in the renormalizable superpotential by the SM 
SU(3) C x SU(2) L x U(l) y vertical symmetry, the U(2) H 
horizontal symmetry and CP conservation are the third 
generation ones and the so called /i-term, 

\tQzU 3 H u + \ b Q 3 V 3 H d + X T C 3 £ 3 n d + fiH u n d . (9) 

We note that, in principle, two other couplings, could 
be allowed in the superpotential: C 3 H U and Q 3 C 3 T> 3 . 
There are different ways to remove this unwanted cou- 
plings. They could be forbidden imposing i?-parity con- 
servation defined as R = (— ) 3 - B + L + 2S ; where B is the 
barionic number, L the leptonic number and S the spin. 
A third possibility would be to extend the \](2) H sym- 
metry to the maximal U(3) ff horizontal symmetry. The 
breaking of the U(3) ff symmetry in the direction of the 
third generation would leave us with our \J(2) H symme- 
try, in such a case this bilinear interaction would not 
be allowed by the U(3) ff symmetry. We also note that 
the couplings in the renormalizable superpotential can- 
not carry complex phases since CP is an exact symmetry 
at this level. Therefore, at tree level the Yukawa matrices 
are generically of the form, 



y 



000 
000 

y 



(10) 



Additionally, trilinear soft supersymmetry breaking 
terms are generated by operators generically of the form, 



E — 



(11) 



— J d 2 9 (T a ^ R + 4> L T a ^) H a , (12) 

where Mp is the flavor breaking scale, o = 1,2 are fla- 
vor indices, 7i a ( a — u,d) represents any of the Higgs 
superfields. Soft supersymmetry breaking mass matri- 
ces can receive diagonal flavor conserving contributions 
of the form, 



1 



d 4 9{ Z^Z' + + ^(t>). (13) 



Additionally non-diagonal flavor violating contributions 
arise from operators generically of the form, 



1 

Ml 



( E 

Z,Z'=S,A 



fa -r-f -r-6,- 
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where a ^= b. Flavor violating supersymmetry-breaking 
fields cannot generate masses for the gauginos. Therefore 
we need to introduce at least one flavor-singlet chiral 
superfield, Q, whose F-term component gets a non-zero 
vev giving masses to gauginos from operators of the form, 



1 

M 



d 2 eg\\ 



(15) 



We note that (Q) breaks supersymmetry but not the fla- 
vor symmetry. We will identify M with the usual su- 
persymmetry breaking messenger scale. We note that 
the messenger scale is in general different from the flavor 
breaking scale even tough the flavor breaking fields are 
supersymmetry breaking fields themselves. The gaugino 
mass generated is given by m~ = (Q) /M. The flavor sin- 
glet superfield responsible for generating gaugino masses, 
G, will couple to matter fields too generating soft trilin- 
ears couplings, 



M I ^ 



(16) 



where k is a real dimensionless coupling determined by 
some unknown underlying renormalizable theory. For 
practical purposes we will assume that k can take ar- 
bitrary values. Soft masses would also be generated by 
operators generically of the form, 



v 



MM F 



(17) 
(18) 



Here flavor indices have been omitted, rj and rj' are also 
real dimensionless couplings determined by the unknown 
underlying renormalizable theory. Regarding the possi- 
ble appearance of D-terms in the scalar potential. D- 
terms would appear when a local symmetry is sponta- 
neously broken by scalar fields, which is not the case for 
the model under consideration. After the XJ(2) H flavor 
and the CP symmetry are broken spontaneously by the 
supresymmetry breaking fields defined in Eq. [7] the fol- 
lowing boundary conditions for the soft trilinear matrices 
are generated at the scale Mp, 



A = A 



crAV*' 

-o-A 2 e^° crAe l * £ 



2o-Ae^ 
1 



(19) 



Km— and the dimensionless parameter a is 



where A 

defined by a — rhj A. The mass parameter rh, defined in 
Eq. [SJ is a new mass scale introduced in the problem by 
the flavor violating susy-breaking fields. We note that 
only one combination of the complex phases in Eq. 1191 
will be transmitted to the Yukawa matrices. For conve- 
nience one can remove some of them from the soft tri- 
linear matrix through a redefinition of the phases of the 
matter fields, even though they will appear in the soft 
mass matrices. Without any loss of generality we can 



adopt a flavor basis where the soft trilinear matrix takes 
the following form, 



A = A 





-a\ 2 



aX 2 aX 2 e~ l ' i 
crA 2ctA 



crA 2 e-^ 2crA 



-i<f> 



(20) 



The phases 7 and <fi are related with the phases in Eq. ^3 
by 7 = ~{<t>f + ®S ~ <i>a- $f) and = (2$ F - $ s ). 
After the U(2) ff flavor breaking flavor violating soft mass 
matrices are also generated. In the flavor basis adopted 
in Eq. 1201 thev take the following form, 



M 



L,R 



mix 



1 + 5pA 2 
pA 3 (2e^ - 
p'X 2 e 1 ^'^ 



) 1 + 5pA 2 

2p'Ae"^ 



') P'A 2 



-i(~f-<t>) 



(21) 
p and 

p' are dimensionless parameters defined by p = in 2 /fh 



Ip'Xe^ 
1 + 5pA 2 

b a - 2$ F ). 



where, m 2 — r/m- and <j)' — 

" •' •* ; ' '- ' - ::1 ^s 

and p' = (r)' /rj)(m/m^). We note that in this scenario 
the amount of flavor violation in the soft mass matrices is 
determined not only by the powers of A in the off-diagonal 
entries but also by the parameters a, p and p' . We note 
that in the limit m — * all the flavor violation will be 
suppressed. There are other interesting limits: if 77' — »• 
the mixing between third and first or second generation 
in the soft mass matrices is suppressed, if rh <ti rhf the 
flavor mixing between first and second generation in the 
soft mass matrices is also suppressed and the sfermions 
masses will be nearly universal, in the case A ~ rh -C rhf 
only the contributions from soft masses to flavor violating 
processes would be suppressed while the flavor violation 
in the soft trilinear matrices could be sizeable. 



III. RADIATIVE GENERATION OF YUKAWA 
COUPLINGS 

In the presence of flavor violation in the soft sector, the 
left and right handed components of the sfermions mix. 
For instance, in the gauge basis the 6x6 down-type 
squarks mass matrix is given by, 



M 2 = 



M 2 DL 



v 2 c 2 yly D (a^ C/3 



(A D cp - fiy^spjv 



+ Vc 



/J,yDS/3)v 

y D yh 



(22) 

where ■M 2 j R and M. 2 Dl are the 3x3 right handed and 
left handed soft mass matrices (including D-terms), 
is the 3x3 soft trilinear matrix, y^, is the 3x3 tree- 
level Yukawa matrix, tan [3 is the ratio of Higgs expecta- 
tion values in the MSSM, p is the so-called mu-term and 
v = swmw/y/2Tra e = 174.5 GeV. M 2 D is diagonalized 
by a 6 x 6 unitary matrix, Z D . The presence of flavor 
violating entries at tree level in the soft supresymmetry 
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breaking matrices will generate one loop contributions 
to the Yukawa matrices. In general, the dominant fi- 
nite one-loop contribution to the 3x3 down-type quark 
Yukawa matrix including CP phases |36j | is given by the 
gluino-squark loop, 



37T 9 ^ 



ac Z (b+3)c B o(mr 



where d c (c = 1, • • •, 6) are mass eigenstates and m~ is the 
gluino mass. Bq is a known function that can be found 
elsewhere in the literature. The radiatively corrected 3x3 
down-type quark mass matrix is given by, 



m D = vc p {y D + y r D ad ). 



(24) 



We note that the effective supersymmetric model pro- 
posed generates an approximately degenerate squark 
spectra. In the squark degenerate case one obtains a 
simple expression for yfj , 



-,;rad _ 2«£ *, . 



/j,y D tan P)F(mr, mr,m~), (25) 



where the function F(x, y, z) is a form factor of the par- 
ticles in the loop with units of [Mass] -2 . F(x,y,z) is 
defined in Ea. ll48l of the appendix. For the soft-trilinear 
texture in Eq. 1191 predicted by our model one obtains a 
simple expression for the radiatively corrected down-type 
quark mass matrix, 



m D = m b 



loY uYe- 1 
—uj\ 2 loX 2uoX 
wA 2 e" i7 2loX 1 



(26) 



Refs. |24ll25ll37| . In this subsection we briefly summarize 
results included in those references. Although not diag- 
onal in the gauge basis, the matrix txid can be brought 
to diagonal form in the mass basis by a biunitary diago- 
nalization, (V^ixidV^. = (m^, to s , m{,). The down-type 
quark mass matrix given by Eq. 1261 makes the following 
predictions for the quark mass ratios to leading order, 



m d 



Ills 



A 2 + 0(A 4 ), 



m s 
mb 



uj\ + 0(X 4 ). 



(29) 



We can relate A and lo with dimensionless fermion mass 
ratios. To first order, 



m d 



1/2 



1/2 



(30) 



Using these relations and the runnin g qu ark masses de- 
termined from experiment, see Ref. |37j | for details, we 
can determine A and lo, they are given by A = 0.209 ± 
0.019 and ui = 0.109 ±0.030. We observe that constraints 
on the supersymmetric spectra can be derived from the 
parameter lo. To assess the viability of the model we 
must check if it is possible for lo to reach the values re- 
quired by the observed quark masses. From Eq. 1371 we 
obtain the following unequality for m~ < rrt^, 



(31) 



Using the measured values v = 174.5 GeV, m^m^MS 
4.2 ± 0.2 and a s = 0.117 we obtain, 





where lo encodes the dependence on the supersymmetric 

spectra. For the case m~r > m~, we obtain 
^ b — g 1 



(27) 



We emphasize that m is not any squark mass scale but 
a new mass scale introduced in the problem by the vevs 
of the flavor violating susy breaking fields, see Eas. 15171 
The parameter rh^ defined as, 

fh b = vc p [y b + LObie'^ - -^-y b , (28) 

is approximately the running bottom mass. Lo b — 
Lofhbl (vcp). The phase e -1 ^ is an overall phase absorbed 
in the definition of fhb in Eq. 1281 which has no observable 
implications. 



LO < l.5C/3 — 



m~ 

/ 



(32) 



Therefore the values of lo required by the measured quark 
masses can be easily reached without any ad-hoc tuning 
in the supersymmetric parameter space. Let us examine 



with some detail the case 
tan (3 = 50, we obtain, 



For large tan (3, 



W «0.03 l™ 



which would require m « ?>m~. On the other hand in the 
opposite gluino mass limit, m~ > 2m~, one obtains, 




A. Quark masses 



For the large tan/3 case , tan/3 = 50, we obtain similar 
results, 



The implications for fermion masses arising from a 
matrix similar to the one in Eq. 1261 were studied in 



lo w 0.03 



In 
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If to ~ m~ we would need to— ~ 20to~. This would imply 

~ 9 9 fe 

m ~ 20m~. The mass matrix in Eq. 1261 is very succesful 

in reproducing the down-type quark mass ratios, but it 
cannot explain correctly the measured mass ratios in the 
up-type quark sector. We will propose a simple solution. 
Let us assume that down and up type quark fields trans- 
form differently under certain Z 2 symmetry. If this was 
the case the fields S, A and T would not generate any 
mixing in the up type quark sector. Let us assume that 
there is an extra C/(2)h symmetric tensor, S' , that gets 
a vev of the form, 



<S' >-- 



X 6 






A 2 



(34) 



If the couplings of S' with the up-type quark fields are 
allowed by the Z 2 symmetry they would induce a soft 
trilincar matrix of the form, 



aX 6 
a\ 2 
1 



(35) 



We note that we do not show any phases in the ma- 
trix A[/. Possible phases in the diagonal entries of Ajj 
are not physical to leading order since they can be ab- 
sorbed through a redefinition of the phases of the matter 
fields. Masses for the up and charm quarks are gener- 
ated radiatively. One can perform an analysis symilar to 
the analysis in the down-type quark sector and obtain a 
simple expression for the radiatively corrected up-type 
quark mass matrix, 



m/ ■ = mt 



ujX 6 
ujX 2 
1 



(36) 



where u) for m~ > to— is given in this case by, 



- * ( — ^ J I — J ( z- 

\mt 07T J \ TO~ J \ TO~ 



and fht is the normalized top quark mass given by, 



(37) 



fh t = vs/3 \^yt + w t (l - -j-yt cot f3) j , (38) 

with cut = ujfht j '(v s p) . For the soft trilinear texture under 
consideration in Eq. [21)1 one obtains the following predic- 
tions for the up-type quark mass ratios, 

^ = A 4 + G(A 6 ), ^=ujX 2 + 0(X 5 ), (39) 
m c m t 

Again we can relate A and uj with dimensionless quark 
mass ratios, to first order, 



m , 



1/4 



m 2 m u 



1/2 



(40) 



Using the run ning quark masses determined from experi- 
ment, see Ref. [33 for details, we obtain A = 0.225±0.015 
and w = 0.071 ± 0.018. We note the similarity in the val- 
ues for A and u> calculated in the up and down type quark 
sectors from Eas. 1301 and 1401 Let us examine with some 
detail the case m~ w m~. Using the measured top quark 
mass, m t = 178 GeV, we obtain the condition, 



lo m 0.02s f3 



TO 



For instance, for large tan /3, tan (3 — 50 we obtain 



u> 



0.02 



This constraint is compatible with the analogous con- 
straint arising from the down-type quark sector, which 
was uj ~ 0.03 (fh/m~), see Eq. I32fl . Therefore no im- 
portant splitting between the sbottom and stop quark 
masses is required for the viability of the model. In- 
deed both constraints could be satisfied simultaneously 
for m~ ~ 1.5to~ which is a non-trivial consistency check 
both of the model and the ad-hoc vevs introduced in 
Eas. 15171 



Radiatively generated CP phase and CKM 
elements 



Finally, one can calculate the CKM mixing matrix. 
This is defined by V C km = V£ f Vf . We have seen in the 
previous section that in the simple model here proposed 
the up-type quark mass matrix is diagonal. Therefore 
the CKM matrix is given by Vckm = Vf,- The diago- 
nalization of the down-type quark mass matrix in Eq. 1261 
lead us to the following expression for Vckm to leading 
order in powers of A, 



1 + 2is 7 Xuj - A 2 /2 -A(H 
-A(1 + 4^A) i(A 2 



- cjA(4 + 2is 7 )) 
+ Au 2 X 2 ) - 1 
2loX 



u;X 2 e- l ~> 
2loX 

LoX 2 (2-e^) 2uoX \-2u?X 2 

(41) 

It is easy to check that this CKM matrix is unitary to 
order A 3 , i.e. V ckm Vckm = I + C(A 3 ). We note that 
the model predicts that to leading order \V US \ = X. The 
measured value of \V US \, |^is| cxp = 0.220±0.0026, agrees 
perfectly with the value of A as calculated from quark 
mass ratios in Eqs. |20l an d 0U1 The model also predicts 
that to leading order uj = \V c b \ /2 \ V US \. Using the mea- 
sured value for \V cb \, \V cb \ cxp = 0.0413 ± 0.0015, we ob- 
tain that uj = 0.093 ± 0.005. We note that this value 
of uj is surprisingly consistent with the value calculated 
from quark mass ratios from Eas. 1301 and!40l Finally us- 
ing these values of A and uj we can predict \V u b\ to be 
(V^l = ujX 2 = 0.0045 ± 0.0003, which is consistent with 
the measured value, |K&| exp = 0.00367 ± 0.00047. Using 
again the previous values of A and uj calculated from \ V US \ 
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and \V c b\ we find the following predictions for \V CS \ and 
|F ti ,|, \V CS \ = 0.9795 ±0.0007 and \V tb \ = 0.9991 ± 0.0001, 
which are also in agreement with experiment. It is a 
trivial check to prove that the angle 7 introduced in the 
parametrization of the CKM matrix given in Eq^J coin- 
cides with the standard definition for 7 = </>2, 



7 = Arg 



V ud V 



lib 



v cd v; b 

The angle 4>i is defined as usual by, 



= Arg 



v td v t 



lb 



(42) 



(43) 



We note that we are using the notation (pi and not the 
usual (3 to avoid any confusion with the supersymmetric 
parameter tan /3, the ratio of the Higgs vevs in the MSSM. 
The angle a = <f>^ can be obtained from, (a + <pi +7) = tt. 
Using our parametrization for the CKM matrix we obtain 
to leading order in powers of A a simple relation between 
the angles <p± and 7, 



h = Arg [2 - e*i) 



(44) 



l xp = 23.3° ± 



Using the measured value of the <f>\ phase, 
3.2° (2cr) we predict to leading order the phase 7 to be 
within the range 7thco = 103° ± 13°. This is far from the 
current la global fit value for 7, j& t = 61° ± 12°. To 
find better agreement with the experimental constraints 
on 7 it is crucial to include in Eq. ^]the next to leading 
order corrections to the unitarity of the CKM matrix. 
We find that the only relevant correction which affects 
the prediction for 7 is the correction to the element Vtd, 



NLO 



td 



(2(1 + 5lu\) ~ e^)cuX 2 



(45) 



Including this correction we predict the phase 7 to be 
within the range 7thco = 91° ± 18°, which intersects with 
the current la global fit value for 7. For this range of 
7 \V td \ is predicted to be, \V td \ = 0.0043 ± 0.0005. Fi- 
nally \V ud \ is predicted to be \V ud \ = 0.9765 ± 0.0006. To 
sum up, in this model not only the first and second gen- 
eration fermion masses but also the CKM phase can be 
generated radiatively in prefect agreement with current 
measurements. 



IV. SUPER CKM BASIS 

Overcoming the present experimental constraints on 
supersymmetric contributions to flavor changing and CP 
violating processes is a necessary requirement for the 
consistency of any supersymmetric model |38j . In our 
scenario, as a consequence of the approximate radiative 
alignment between Yukawa and soft trilinear matrices 
there is an extra suppression of the supersymmetric con- 
tributions to some of these processes. Therefore for calcu- 
lational purposes it is convenient to rotate the squarks to 



the so-called superCKM basis where this radiative align- 
ment mechanism is manifest. 

The superCKM basis is the basis where gaugino ver- 
tices are flavor diagonal H^,^3j^|. In this basis, the 
entries in the soft trilinear matrices are directly propor- 
tional to the corresponding contributions to flavor chang- 
ing proceses. For instance, the soft trilinear matrix Ad 
in the superCKM basis is given by, 



a SCKM 



= (V d jA D V d R , 



(46) 



where R are the down type quark diagonalization 
matrices. The soft trilinear matrix Ad is given by 
Eq. 1201 The Yukawa diagonalization matrices are given 
by = Vckm in Eq. |2] while Vj| is completely deter- 
mined to obtain real mass eigenvalues after the diagonal- 
ization of ni£). We obtain, to leading order in A, 



Im[A s D CKM ] = { _ )Ab 
while the real part, Re[Afp KM ], is given by, 



2ws 7 ctA 3 s 7 ctA 2 
2s,/,wA 
2s ct ,uiX Scf, 



2s 7 cto;A 3 
s 7 ctA 2 

, SCKMl 



(47) 



A,, 



aX 3 



4ctwA 3 (c 7 - 1) ctA 2 (c 7 -2) 
+ 1) Act 2X(c l f > uj — a) 

ctA 2 (c 7 + 2) 2A(c0oj — ct) C0 + 8u;aX 2 



4ctwA 3 (c 



(48) 

We note that the entries (2, 1) and (1, 2) contain an ad- 
ditional suppression factor luX compared with the soft 
trilinear matrix in the flavor basis, see Eq.[3Ul This sup- 
pression is a consequence of the radiative alignment be- 
tween Yukawa and soft trilinear matrices. It is convenient 
when calculating supersymmetric contributions to flavor 
violating processes to use the parameters (Sfj)LR defined 
as, 



(4) 



LR 



/-aSKM\ 
VCf3{A D )i2 



(49) 



For consistency we also need to calculate the down-type 
squark soft mass matrices in the SCKM basis. For in- 
stance, for the left-handed soft mass matrix we obtain, 



(M 2 Dl 



\SCKM 



(50) 



and analogously for the right handed soft mass matrix. 
Assuming the soft trilinear texture from Eq.|^we obtain 
for Rc[(M 2 Dl ) sckm ], to leading order in A, 



1 yA 3 y'X 2 
yX 3 1 -2 C< pp'X 
y'X 2 -2c^p'X (1 + (8c p'w - 5ctA 2 )A 2 ) 



Here m~ , y and y' are defined as, 



m 2 r (l + (l + 5p)A 2 ), 



(51) 



(52) 



y = (ip'u^Cj, - 2c0) + p(c 2 0' - 2c 7 )), (53) 



y' = ( c (4>-i) ~ 2c 4>)p'- 



(54) 
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We note that if the gluino mass is of the same order than 
the squark masses, m; ~ fhf, fhf ~ fh and 77' ~ 77 we 
expect that p ~ p'. In that case the coefficient y simplifies 
to y ps p(c20' — 2c 7 ) since u; » 2A 2 . Furthermore if p'/p = 
(77' '/ '77) (to 2 / 'fhm^) the limit p' 3> p would correspond to 
77' ^> 77 or m -C m 2 or to^ -C to 2 , if this is not the case 
we would expect the p term to dominate. The imaginary 
component, \m\{M. 2 D ) SCKM ], is given to leading order 
in A by, 



-z\ 3 -z'\ 2 

-2s p'A 
2s lb p'\ 



(55) 



where, 



z = (Ap'usCjS^ + p(s 2 0' - 2s 7 )), 
2' = (2s0 - s ( 0_ 7 ))p'. 



(56) 
(57) 

If p ~ p' the coefficient z reduces to z w p(s2^' — 2s 7 ) 
since w ~ 2A 2 . Assuming the soft trilinear texture from 
Eq.EIlwe obtain for Re[(Mf 3fl ) SCKM ], to leading order 
in A, 



1 



Here. 



rA 3 
■A 3 1 

'A 2 -2c p'A (1 + 



771^1 + (l + 5p)A 2 ), 



r'A 2 
2c^p'A 

'uj - 5ctA 2 )A 2 ) 



(58) 



(59) 



r = (4p'w(c (7 _0) + 2c^) + p(c 2 0' - 2c 7 )), (60) 
r' = (c (7 _ 0) +2c )p'. (61) 

Again if p ~ p' the coefficient r reduces to r w 
p(c20' — 2c 7 ) since u> « 2A 2 . The imaginary component, 
Im[(A4|j ) KM ], is given to leading order in A by, 



where, 



-a 3 -i'A 2 
t\ 3 -2s p'A 
i'A 2 2s p'A 



£ = (S20' - 2s 7 )p, 
i' = (s( 7 _0) - 2s )p'. 



(62) 



(63) 
(64) 



It is also convenient when calculating supersymmetric 
contributions to flavor violating processes to use the cou- 
plings {5f,) LL defined by, 



1 M 2 \SCKM 
(Oij)LL - 2 ' 



One can define analogously the (SfAjm couplings. 



(65) 



V. SUPPRESSED CONTRIBUTIONS TO EDMS 

In an unconstrained minimal supersymmetric standard 
model (MSSM) the generic contribution to the neutron 
EDM [lH . l26l 23 is around eights orders of magnitude 
larger than the SM contribution, i.e. about four or- 
ders of magnitude above the current experimental con- 
straint ■ This is the so called SUSY CP problem or to 
be more especific the flavor conserving SUSY CP prob- 
lem. The disparity between the current experimental 
constraint and the generic supersymmetric contribution 
in the MSSM is due to the, in principle, allowed presence 
of CP phases in the superpotential and in the soft super- 
symmetry breaking sector. Numerous papers have exam- 
ined this topic in the context of supersymmetric models 
|44j and a few solutions have been proposed, which were 
summarized in the introduction. We will explain with 
some detail how generic supersymmetric models for ra- 
diative mass generation can ameliorate this problem. We 
will analyze separately the one-loop, two-loop and higher 
order contributions to EDMs. 

Interestingly, the one loop supersymmetric contribu- 
tions to EDMs always appear as combinations of six pos- 
sible physical phases of the generic form [lj, |4{J , 



Arg(A*TO~), Arg(B>TO T ), 



(66) 



where A are first generation flavor diagonal trilinear 
soft supersymmctry breaking parameters, to are gaug- 
ino masses, B is the bilinear soft supersymmetry break- 
ing term and pi is the superpotential bilinear term. In 
the special case of universal soft supersymmetry break- 
ing terms these reduce to only two physical phases. First 
let us focus our attention on the term Arg(_B* pm~) . The 
tree level gaugino masses are flavor conserving param- 
eters generated by the supersymmetry breaking flavor 
singlet field Q as we explained before and as a conse- 
quence cannot carry complex phases, which are linked to 
flavor breaking vevs. The p term is allowed in the CP 
conserving superpotential at tree level. This parameter 
is linked to the flavor blind operator TL u 'Hd 1 which obvi- 
ously cannot carry CP phases at tree level. For the same 
reason the bilinear soft supersymmetry breaking term, 
B, is also a real parameter since the term h u hd in the 
scalar potential is also a flavor singlet. Regarding the 
contributions of the form Arg(yl*TO~). In the case of the 
neutron and mercury EDMs the relevant terms arise from 
the up and down quarks EDMs which are Arg(^4*m~) 
and Aig(A^m~). We have seen that in our model the 
3x3 matrix Ay is diagonal, see Eq. 123 The diagonal 
entries carry no CP phases. Even whether they existed 
they could be absorbed in a redefinition of the phases of 
the up-type matter fields. Furthermore, the entry (11) 
of the 3x3 soft trilinear matrix in the down-type squark 
sector corresponding to Ad in Eq. 1471 is in general real. 
Therefore all the one-loop contributions to the EDMs in 
this model are exactly zero. 

Regarding the two loop contributions to EDMs. It 
has been pointed out that the two-loop supersymmetic 
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contributions can also constraint the supersymmetric pa- 
rameter space, even though not so severly as the one-loop 
contributions 

El El 49] . For instance, two- loop contri- 
butions of the Barr-zee type to dd exist with an exchange 
of stops and the CP-odd Higgs. This contribution is of 
the form E3, 



dd 



-a e m d mt ^sin(26> t ) 

777T 3 V 



sin S t C(t u t 2 , A) 



2-ioop "" " s 20 m A 

(67) 

where St = Arg[A t + cot/3/x*], 9 t is the stop mixing, 
v = 175 GeV, A is the CP-odd Higgs and C stands for 
a dimensionless two-loop form factor which can be found 
in Ref . [53 . For the large tan (i case under consideration 
St ~ ArgL4 t ] = since A t and /i are real parameters. 
On the other hand there is an analogous contribution 
from bottom squarks that is proportional to ImL4be I<5i '] 
with 8b = ArgLAb + tan/3^*] For the large tan/3 case un- 
der consideration the second term in 8b would dominate 
in general and we obtain that the contribution is pro- 
portional to — Absmcf), as can be seen from Eq. |27| We 
obtain that for CP-odd Higgs masses above 1 TeV the 
resulting contribution to the neutron EDM is below the 
current experimental constraint, \d n \ < 10 -25 e.cm. 



One may wonder if the previous arguments for the can- 
cellation of the one loop contribution to EDMs could be 
extended in a variant of this model to not only suppress 
but cancel the two-loop contributions. We note that if 
the flavor model generates hermitian soft trilinear matri- 
ces the two-loop contributions, which are proportional 
to the phases of Ab and A t , would be zero. It may 
be possible in principle that small phases are generated 
in the "flavor conserving" parameters in the lagrangian, 
as the /i-term, the B parameter or the gaugino masses. 
Nonetheless, we note that many of the higher order oper- 
ators which could contribute to the radiative corrections 
to the n term have to be flavor conserving operators of 
the form: cx Z^ ab Z ab (Z = S,A) or cx 7* a T a , etc,--. Af- 
ter the breaking of the flavor symmetry these operators 
cannot generate complex phases since the presence of a 
complex phase would be an indication of flavor violation. 
To sum up, because of the intrinsic flavor off-diagonal 
nature of the CP violating phases in this model, the one 
loop contributions to EDMs are zero and two and higher 
order contributions are suppressed below experimental 
limits. 



VI. CONTRIBUTIONS TO DIRECT AND 
INDIRECT CP VIOLATION IN THE KAON 
SYSTEM 



The measure of indirect CP violation in the Kaon sys- 
tem is given by the parameter e defined by, 



A(K L 



iir/4 



A(K< 



V2A 



rriK 



-Im [Mi 



(68) 



where Am a: is the K^Ks mass difference, Ai\ 2 = 
M{K°) = (K° \H* S S=2 \ K°) is the K°K° mixing am- 
plitude and Htff =2 is the effective AS = 2 hamilto- 
nian. The parameters ArriK and |e| have received con- 
siderable attention in supersymmetric models since their 
measured values have been known with good precision 
for long time (H E3, EH- We will separate the SM 
and supersymmetric contributions to the mixing ampli- 
tude in the form Mi 2 = Mf 2 SY + We define 4> e 
and <p' t as the phases of the SM and the supersymmetric 
contributions respectively, Mf 2 SY — l-Mf^^e*^ and 



Mi2~ — IMf^le*^. It is also convenient to introduce 
the ratio R K = \Mf$ SY \/\Mff\. This ratio and the 
complex phase (j)' e are constrained by the experimental 
measurements of Auik and |e|. We can expand in pow- 
ers of Rk to obtain the following expression for the new 
physics contributions to Amx and |e|, 



ArriK — Ai 



,SM 



Am| M 



C SM| 



C SM| 



Rk 



Rk 



cos(2$) 
cos(20 e ) 

sin(2#) 
sin(2<^ e ) 



(69) 
(70) 



We note that ArriK has been measured with an uncer- 
tainly of approximately 0.2%, Am^ p = (3.490±0.006)x 
10~ 12 MeV. In the SM roughly 70% of the measured 
ArriK is described by the real parts of the box diagrams 
with charm quark and top quark exchanges. Some non- 
negligible contribution comes from the box diagrams with 
simultaneous charm and top exchanges while approxi- 
mately the remaining 20% of the measured Avcik is at- 
tributed to long distance contributions. On the other 
hand these are potentially sizeable and up to date incal- 
culable. While a precise prediction is not possible, the 
observation is roughly compatible with the SM expec- 
tation. Assuming that the supersymmetric contribution 
saturates one half of the experimental measurement and 
using Eq. 1691 we can obtain an approximate upper con- 
straint on Rk, 



Rk cos(2cf>' e ) < icos(20 e ). 



(71) 



|e| has been measured with an uncertainity of approxi- 
mately 0.6%. A fit to the K — > tttt data yields |e| = 
(2.284 ± 0.014) x 10~ 3 . The calculation of |e| is also af- 
fected by large distance corrections. It is usual when cal- 
culating |e| to input the experimental measurement for 
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Atok in the denominator of Eq. [58] We will allow again 
for the supersymmetric contribution to saturate one half 
of the experimental measurement. The resulting con- 
straint on Rk and <j)' £ can be expressed in the form, 



R K sm(2cf>' e ) < isin(2#). 



(72) 



Using the usual Wolfenstein parametrization of the CKM 
matrix to second order in powers of A we obtain that 
the dominant SM contribution to M± 2 is proportional 
to (V cs V c * d ) 2 . This is given by V cs V* d = -A(l - £)(1 + 
A 2 X 4 (p — irj)). This can be written as V cs V* d « — A(l — 

^-)e~ l ^ e where 4>e is defined as <p t = tan _1 (?7A 2 A 4 ). This 
is a very small number. Using updated extractions of A, 
r\ and A we obtain <f) e w 0.03°. Therefore the constraints 
from Eas. 1711 and 1721 can be written in the form, 



R K cos(20' e ) < 0.5, 
E^sin^) < 6 x 10" 4 . 



(73) 
(74) 



Next we need to calculate the expressions for Rk and <f/ e 
in our model. We will see that <p' e in the model under con- 
sideration is completely determined given the measured 
value of 7 and the quark mass ratios. Therefore the previ- 
ous constraints in Ea . 1731 and 1741 will translate into lower 
bounds on the squark mass spectra. The supersymmetric 
contribution to M \ 2 contains contributions proportional 
to the different 5 couplings in the soft supersymmetry 
breaking sector. There are qualitatively four different 
contributions to A4f 2 SY that in our model, to leading 
order in powers of A, take the following values, 
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+ (*Hi)ia = 28¥ K (-c 2 ~ l + --i S2l ), (75) 
5 5 

(<Wl2(<Wl2 = ^r(^ c 27 ~ 4 ~ * S 2 7 ), (76) 



where, 



^ 2 R = 8 W W 



2 — 2 

b b 



and, 



(St 



LL)12 "T {&RR)12 
d \ /cd 



26}i(e- 2lrt >' -2e"^) 2 , 



(<*ZI,)l2(<W) 



12 



J LL\ 



— Ot 



2e" 17 ) 2 , 



where, 



511 = P 2 A 6 



(77) 

(78) 
(79) 

(80) 



For instance, we noted in Sec. IIII Xj that for large tan/? 
and to— ~ m~ the parameter to (new mass scale intro- 
duced by the flavor violating susy breaking fields) is re- 
quired to be approximately to sa 3m^. Therefore the 
LR couplings in our model have an important additional 
suppression factor, S 1 LR jb\ L oc c|w 2 u 2 /m~. As a conse- 
quence in this model the Sll and Srr couplings dominate 



the contribution to A4f 2 SY . All the Sll and 5rr contri- 
butions to the Wilson coefficients can be added up in a 
simple expression. We obtain, 



^SUSY 



a 2 s m K f K X K 



l 12 



60to2- 

b 



xfl(i) 



(sL)Ms d RR ) 12 c(x)] 



(81) 



Here X K is a dimensionless factor defined as X K = 
m Kl( m s + m d) 2 , numerically Xk ~ 4.07, fx is the K- 
meson decay constant and x is the gluino-squark mass 
ratio squared, x = m^/mA. C(x) and D(x) are given by 

C{x) = C f f{x) + C g g{x) and D(x) = Dh{x) where f(x), 
g(x) and h[x) are dimensionless form factors defined in 
Eqs. 11491 11501 and 11511 of the appendix. The functions 
f(x), g(x) and h(x) have been conveniently normalized 

so that in the limit to— ss m~ they tend to 1. We note 

g b •> 

that in our approach the constant coefficients Cf,C g and 
D absorb the dependency on the method of calculation 
of the hadronic matrix elements as well as the renor- 
malization effects on the Wilson coefficients. Following 
Ref. [5l|, where lattice QCD methods were used to calcu- 
late the relevant hadronic matrix elements and including 
NLO renormalization effects to the Wilson coefficients, 
we obtain the following numerical values, Cf rs 9.13, 
C g fa 0.75 and D m 0.002. We note that the naive 
vacuum insertion approximation at tree level gives the 
values C f « 1.8, C g w 0.067 and D w 0.0055, which are 
significantly different. We also note that these constant 
parameters, Cf,C g and D, do not depend on the flavor 
mixing structure in the soft supersymmetric breaking sec- 
tor. From the numerical values of these coefficients we 
note that the contribution of the form SllSrr in Ea. 1791 
dominates the supersymmetric contribution. Finally, us- 
ing the well known expression for the SM contribution to 
the amplitude, see for instance Eq.(3.39) in Ref. |^3, we 
can write the ratio of the supersymmetric contribution 
over the SM amplitude in the form, 



R K ^2^f(x)p 2 \ e £-cos( 1 



b 



1/2 



where, 



a s irX K C 



1/2 
7 



666 TeV. 



(82) 



(83) 



Here Bk = 0.85 ± 0.15 is a renormalization group invari- 
ant form of the B parameter arising from the hadronic 
matrix element, m c is the charm quark mass and 771 is a 
short distance QCD correction factor, at NLO 771 is given 
by Tji — 1.38±0.20 . Using for A the value determined 
from CKM elements, A ~ 0.22, and for 7 the la global 
fit 7fit = 61° ± 11° we obtain, 



Rk 



10 TeV 



b j 



h(x), 



(84) 
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The phase of the supersymmctric contribution, <j)' tl can 
be calculated from Eq. This phase depends strongly 
on the phase <fr', which is not constrained by the CKM 
matrix. For instance, for cf>' = we obtain, 



tanf 



C2 7 ) 



(85) 



For 7 = 60° we would obtain cj)' e = 0. If this was the case 
the measurement of |e| would not constraint the super- 
symmetric spectra, see Eq. 1741 The only phase indepen- 
dent constraint comes from the ArriK measurement, see 
Eq. [73] Let us analyze with more detail the large tan (3 
case. For tan (5 = 50 and m~ fts we noted in Sec. IIII Al 
that for the hierarchy of flavor breaking vevs postulated 
in Eas. 15171 we need to have rh ss 3to^, i.e. p w 9. There- 
fore using Eas. 1731741 andl%2lwe obtain for tan/3 = 50 the 
constraints, 



> 130 TeV, 
rrt£ > 3670 TeV, 



Am K , 

\e\ (^0). 



(86) 
(87) 



To obtain the second constraint we assume that the phase 
<f>' is arbitrary. If this phase was zero as we mentioned 
this second constraint would not be effective. We would 
like to emphasize that these constraints are non-generic. 
They only apply for the case large tan fj case study and 
the particular texture considered in this paper. This case 
study must be considered a worst case scenario for this 
kind of models. Several variants of this model may allow 
us to lower considerably these bounds. For instance we 
could lower the bounds at the price of less predictivity in 
the fermion sector by increasing the number of parame- 
ters in the flavor breaking sector. 



The measure of direct CP violation in the Kaon system 
is given by the parameter e'/e. Direct CP violation orig- 
inates from direct transitions of the CP-odd state into 
the CP-even irir final state. The direct CP violation in 
the neutral K — ► tttt decays can be described through the 
ratio, 



- - e 1 * "' 
e 



V2\ 



ImA 2 ImAp 



ReA 2 ReA 



where Aq. 2 are the isospin amplitudes for the AI — 
1/2, 3/2 transitions. lmAo t2 are calculated from the gen- 
eral low energy effective Hamiltonian for AS — 1 tran- 
sitions jjsjl . w = ReA 2 /ReA and $ is a strong phase 
shift difference between the two amplitudes. In 1999 the 
NA48 experiment H3 at CERN and the KTeV H| ex- 
periment at FNAL demonstrated that this observable is 
actually different from zero as expected in the SM. The 
present world average is [5(j, 



Re 



(16.6 ± 1.6) x 10" 



(89) 



There is no simple approximate expression for the SM 
contribution to e'/e. This calculation is affected by large 
hadronic uncertainities. It has been recently pointed out 
[57| that to lowest order (in 1/N C and in the chiral ex- 
pansion) Re [e'/e] is governed by the competition between 
two different decay topologies and suffers from a strong 
cancellation between them. Nevertheless to higher orders 
chiral loops generate an enhancement of the isoscalar am- 
plitude and a reduction of A 2 . Taking this into account, 
and following Ref. [57|], the latest SM prediction is: 



Re 



(19 



SM 



+17 

-18 



x 10" 



(90) 



Therefore even though the SM prediction is consistent 
with the measurement, it does not allow us at present 
to perform stringent tests of the CKM mechanism of CP 
violation. For a 2003 review of several calculations see 
the Refs. H|. 

Since the 1999 measurements, the parameter e'/e has 
received considerable attention in the context of super- 
symmetric theories [59^ . The dominant supersymmet- 
ric contributions to e'/e come from the chromomagnetic 
operators like O g , O g = g s /(l6TT 2 )d L cr^t A s R G^. The 
Wilson coefficient C g corresponding to this operator is 
given by, 



c g = (-) 



2m~ 



4to, 



N(x)(S[ 



LL)l2 



M(x)(6t R ) 12 



(91) 

where N(x) and M(x) are dimensionless form factors de- 
fined in Eas . 1 1 54l and l 1 5 5l of the appendix. Taking into ac- 
count that the relevant hadronic matrix element is given 
by, 



((nn) I =o\O g \K°) = 



3 11 



2 16tt 2 F3 



(ll) 2 n 



(92) 



with Ftt — 131 MeV and where the Bq factor is not well 
known, Bq = 1 — 4 we obtain that the total LR super- 
symmetric contribution to Re [e'/e], can be conveniently 
written as, 



Re 



where 



LR+RL 



11 V3 



-) \lm[(6 d LR hi-(S d LR y i2 ]\N(x), 



rC IT 



64tt |e|ReA° F T (m s + m d ) 



(93) 



a s (mr)T]B G (94) 



We used for w and ReAo experimental values w « 1/22 
and ReA = 3.326 x 10 -4 MeV. For the rest of param- 
eters we used ttik = 490 MeV, — 140 MeV. r\ is a 
well known dimensionless strong coupling renormaliza- 
tion factor defined in Ref. [6(j. Taking into account the 
important uncertainities in the current determination of 
the Bq factor and the lighter quark masses we obtain 
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the estimate = 100 - 600 TeV. We note that for the 
model under consideration 



Im[(^ R ) 12 ] = Im[(^ R ) 21 ] = {-)2uj S ,cpX- 



v m 



m~m~ 

b b 



For the worst case scenario, r\ t i = 600 TeV, and assum- 
ing that the supersymmetric contribution saturates the 
experimental measurement we obtain the constraint, 




N(x) < Re 



(95) 



|6lL l62| , that the CP symmetry is significantly violated 
in the B sector, in agreement with Standard Model pre- 
dictions, providing a confirmation of the so called CKM 
paradigm [63, |64j . This fact does not rule out the possi- 
bility that the effects of CP phases of a different origin, as 
for instance the phases in the soft supersymmetry break- 
ing sector |6jj |66J, |67j , could manifest in the near future 
|68j through other CP violating observables, especially 
penguin dominated amplitudes such as B — ► <f>K° ,r)'K s . 
In this section we will study the constraints that the 
currents measurements of CP asymmetries in several B- 
decays imposse on the model under consideration. 



m~, exam- 

d 1 



For the large tan j3 case, tan (3 — 50 with m~ 
ined in Sec. IIII Al the hierarchy fh 

In this case using for A and u> the values determined from 
the CKM elements and for 7, 7 — 60° we obtain, 



3m~ was required. 



m~> 1.3 TeV (e'/e). 



(96) 



Let us analyze separately the size of the LL and RR con- 
tributions. We see that because of an extra m s /m^ sup- 
pression factor the LL+RR contribution to e' is much 
smaller than the LR contribution. If the gluino mass is 
of the same order than the squark masses (which is re- 
quired to maximize the loop generated quark masses), 
and p — p 1 we obtain a simple approximate expression 
for the LL and RR couplings, 



Im[(<5 LL ) 12 ] - Im[(5 RR ) 21 ] = (2s 7 - s 2 ^)X 



m 



The total LL+RR supersymmetric contribution to 
Re [e'/e], can be conveniently written as, 



Re 



Im[(5 d RR hi (6 d LL y i2 ]M(x) 



LL+RR 



(97) 

Using for A the value determined from the CKM ele- 
ments and assuming that the supersymmetric contribu- 
tion saturates the experimental measurement, we obtain 
the constraint, 



(2s, 




< Re 



(98) 



cxp 



For the large tan f3 case, tan f3 = 50, with m~ ~ m~ and 
m & 3mr we obtain rrtr > 220 GeV. 

b b 



VII. CONTRIBUTIONS TO CP ASYMMETRIES 
IN THE B SYSTEM 

The CP violation measured in neutral K meson decays, 
taking into account current experimental and theoretical 
uncertainities, can be simply explained with the CKM 
phase. B factories have verified, especially through mea- 
surements of the CP asymmetry in the Bd — > ipKs decay 



A. CP asymmetry in B — > tpKs 

The B°B° mixing amplitude is defined by the matrix 
element of the effective AB = 2 hamiltonian as Mb = 
Mi2(B°) = (Bd \H§P= 2 \ B d ). The phase of the mixing 
amplitude is related with the mixing CP asymmetry in 
the decay B — > tpK s by, 



sin [Arg [Mi 



(99) 



According to the most recent averaged experimental re- 
sults of Babar and Belle S i>Ka = 0.736 ± 0.049. This can 
be simply accounted to date with the CKM phase. If fu- 
ture measurements reduce considerably the experimental 
uncertainity in S^k s i there is hope that deviations from 
the SM prediction could be elucidated. It is convenient 
to separate the SM and supersymmetric contributions to 
the mixing amplitude in the form Mb — Mf M + Mf USY . 
We also find convenient to define <pi and <f>i as the 
phases of the SM and supersymmetric component of the 
amplitude respectively, i.e. Mf M = e l2( ^ 1 \M 



f M \ and 



M 



SUSY _ 



3 i2^|^SUSY|_ Since the gM p rec jiction can 



account perfectly for the experimental result we expect 
that the supersymmetric contribution is a small correc- 
tion and expand the expression for the CP asymmetry in 
powers of the ratio = \Mf VSY \/\Mf M \. To leading 
order, 

S i ,K 3 = sin2^i[l - sin(2(0! - + erin(2#).fy. 

( 10 °) 

It is known that in the absence of new physics contri- 
butions the SM CP phase can account for the present 
experimental results for S^k s ■ From Eq. 11001 we can ob- 
tain constraints on R^ and <f>\. Assuming that the new 
physics contribution saturates a 50% of the experimental 
uncertainity we obtain to leading order in R^, 



R^ sin(2(^ - <£i)) < 0.5 



(101) 



We note that even in the limit where the complex phase 
of the SUSY amplitude goes to zero the mixing CP 
asymmetry, S^,k s ? is affected by the SUSY contribu- 
tions through their effects on the absolute value of the 
amplitude, S^ Ks = sin20 1 M [l - sin(20f M )i?^]. The 
mass difference in the BdB d system, Amd — 2Abs[Mb], 
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also puts a stringent constraint on R^,. Amj is an 
observable well known experimentally, to the level of 
1.5%. The experimental measurement yields Am d = 
(3.22 ± 0.05) x 10~ 10 MeV. The SM prediction for Am d 
is about Amf = (2.9 ±2.2) x lO" 10 MeV. We note that 
even tough the theoretical uncertainity is about 75% the 
central value is only 10% from the central experimental 
value. Assuming that the supersymmetric contribution 
saturates a 50% of the experimental measurement and 
expanding Mb in powers of we obtain the constraint, 



ify cos(2(# - 4>i)) < 0.5, 



(102) 



We note that we have assumed that (</>i — </>[) ^ ±7r/4. 
If this was not the case, then the second order term in 
the expansion would be dominant and we would obtain 
a milder constraint, i? 2 < 0.5. 

Next we need to calculate the expressions for R^, and 
<j>i in our model. The previous constraints in Eg. 110 H and 
11021 will translate into lower bounds on the squark mass 
spectra. Again, we identify four qualitatively different 
contributions to A^f USY , that in our model, working to 
leading order in powers of A, take the following values, 

(*L)i3 + (4J?3 = 2<Si 3 R (3 + 2c 2 - 2is 2l ) (103) 
(S d LR )i2(S d RL ) 12 = SH(2d-5-is 2l ) (104) 



where, 



rl3 _ , 4 2 v m 
b b 



(105) 



and, 



(Sll)i3 + (Srr)i3 = 2# 3 L e 2i ^>(l + 4e 2 *),(106) 



13 2i 
LL e 



where, 



A 13 



p X . 



-^(l-4e 2 ^), (107) 



(108) 



We see that in this model the LL and RR delta couplings 
in general are expected to dominate the contribution to 
Im[A1b USY ] since the LR couplings contain an additional 
suppression factor, CpV 2 /m~. It is possible to add all 

the dominant Sll and S RR contributions to the Wilson 
coefficients to give a simple approximate expression for 
the supersymmetric contribution to Mb- We obtain, 



M SVSY 



a 2 s m B f 2 B X 2 B 



60m 2 - 

b 



xD(x) 



(109) 



Here X B is a dimensionless factor defined as X B = 
nig / (mb(mb) + mdirrib)) 2 , numerically X B ~ 1-08. x 
is defined as x — m~/rr£. C(x) and D(x) are dimen- 
sionless form factors which were already introduced pre- 
viously in Sec I VI Al as C(x) — Cff(x) + C g g(x) and 



D(x) — Dh(x). f(x), g(x) and h(x) arc dimension- 
less form factors defined in Eqs. 11491 11501 and 11511 of 
the appendix. The constant coefficients Cf, C g and D, 
like in the AS = 2 case, absorb the dependency on the 
method of calculation of the hadronic matrix elements 
as well as the renormalization effects on the Wilson co- 
efficients. We have evaluated Cf, C g and D following 
Ref. [6^| where lattice QCD methods were used to cal- 
culate the relevant hadronic matrix elements. We obtain 
that D « C 9 « C/ and C f « 7.33. Had we used the 
naive vacuum insertion approximation at tree level we 
would obtain, C f = (48V| + 9)/(27X|) w 2.06 which 
would imply an underestimation of the dominant term 
by a factor of order 1/10. The coefficients Cf,C g and 
D also depend on the scale of the supersymmetric spec- 
tra. We have calculated our numerical values at the scale 
M$ = 1 TeV using the renormalization factors given in 
Ref. The contribution of the form 8ll$ RR in Ea. 11081 
clearly dominates the supersymmetric contribution. We 
use the well known expression for the SM contribution to 
the amplitude Mb, see for instance Eq.(3.60) in Ref. [52| . 



r 2 

M 8M = F 



12tt 



VBB B jlm Bd m 2 w \V td V; b \ 2 S(x t ), (110) 



where r/ B — 0.55(1) is a QCD correction factor, B Bd 
is a renormalization g roup invariant parameter avail- 
able in the literature [7C| . B Bd = 1.30 ± 0.12, and 
S(xt) is a dimensionless form factor given by S(xt) = 
2.46(m t /170 GeV) 152 . We can write the ratio of the su- 
persymmetric contribution over the SM amplitude in the 
form, 



7? I ~ ^ 

I LL+RR ~ ^ 

b 



P ' 2 A 4 /W(17-8c 27 ) 1 / 2 , (111) 



where, 



a s wX B C 



1/2 

/ 



V5G F \V td V tb \ B B ,2 VB ,2 m w S^ 2 {x t ) 



78 TeV. 



(112) 

We will substitute in the expression II 111 the value of A 
determined from the quark data and for 7 we will use the 
central value of the la global fit. 7 = 60°. For the large 
tan P case, tan/3 = 50, examined in Sec. 1111 K\ assuming 
that p' ~ p rs 9 and m~ ss m~ the constraints in Eq. 11011 
and 11021 reduce to, 




£2(^-^1) < 0.5, 



S2(^-0O 



< 0.13, 



(113) 
(114) 



The phase 4>\ can be calculated from Ea. 11071 We obtain, 

4S2</> + S2( 7 _0) 



tan(20i) = 



4c 



2<p ~ C2( 7 - 



(115) 
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The phase <f> is not constrained by the quark masses and 
mixings. The value of (/)[ depends strongly on the phase 
<f>. For <f> = and 7 = 60° we obtain </)[ = 10.9° while 
for tf> = 30° and 7 = 60° we obtain (j>[ = 70.9°. When 
cos(2(<^ — <pi)) w 1 the strongest phase independent con- 
straint on the squark spectra comes from Eq. 11131 







D 



> 



> 



1000 TeV {Am d $ x » <^), (116) 
1000 TeV (S^K a <t>\ ~4>i± tt/2). (117) 



On the other hand if // -C 1 these constraints on the 
squark spectra would be milder. This will happen for 
instance when 77' <C r/. 



B. CP asymmetry in B — > c/>A's 



between the experimental measurement and the SM pre- 
diction, i.e., 



2 s<l}NP c 5 C2^R^ < 0.40 ±0.26 



(119) 



where for S24 >1 we used the experimental value of S^Ks > 
S^,Ks — 0.736 ± 0.049. There are two basic contribu- 
tions to the supersymmetric amplitude: the contribu- 
tions coming from the 5lr couplings and the contri- 
butions coming from the Sll and Srr couplings. We 
will first obtain a simple expression for the <5lr con- 
tributions. Wc note that at the susy scale there is 
only one Wilson coefficient which contains the coupling 
{S'[ R )23, the chromomagnetic operators like O g , O g = 
g s / {\&TT 2 )d L a^t A b R G A v . The Wilson coefficient C g cor- 
responding to this operator at the susy scale is given by, 



The latest results from BELLE collaboration for the 
time dependent CP asymmetry coefficient S^k derived 
from the combined cj)K° dataset are [7l| S^§ LLE = 
+0.06 ± 0.42 while the latest results from BABAR col- 
laboration for the same coefficient are [Z3 S$k BAR = 
+0.50 ± 0.32. Combining the results from both experi- 
ments one obtains the world average, cjBABar+belle _ 
+0.34+0.21 Taking into account that the SM pre- 
diction for the time dependent CP asymmetry is Sqk = 
sin(2i/>i) = 0.726 ± 0.037 (here we used the world aver- 
aged CP asymmetry determined from charmonium final 
states) the current world average seems to differ from the 
SM expectation by about 2a level. Therefore this pro- 
cess is one of the best candidates for the manifestation 
of new physics in the quark sector. Sqk in the context 
of supersymmetric theories has received considerable at- 
tention recentiypl EE EE E3, ES EE E3, EE IzS E3, E3, 

EllIEIMIIEEismce any clear deviation from the SM 
prediction would imply the existence of new CP phases 
other than the CKM phase. 

The total decay amplitude, A^Ks =< - 



4>K s \H c ll =1 \B° > can be written in the form, 
A 4>Ks = Af£ s + Af^ s /. Additionally one can 
parametrise the SM and supersymmetric contributions to 
the amplitude in the form, Af^ Y = |A^ SY |e^ NP e i5suSY 



and Af£ s = \A s ^ s \e l5sM . where N p is the CKM-likc 
complex phase of the supersymmetric contribution 
and #sm and ^susy are the SM and supersymmetric 
CP conserving strong phases respectively. Assuming 
that the susy contribution to the amplitude is smaller 
than the SM one, and expanding in powers of the 
ratio R<f, = |^0Xs Y l/l^™ s l it; is possible to obtain 
approximate expressions for the direct and mixing CP 
asymmetries |74j, 187] . To leading order in R4, , 



S<j>K s — s 2</> 1 + 2s l p Np CsC2 l p 1 R4 



(118) 



where <5 is the difference of supersymmetric and SM CP 
conserving strong phases, 6 — (Ssm — ^SUSy)- We will 
constraint the supersymmetric contribution to S<f,K s as ~ 
suming that this contribution accounts for the difference 



C 9 (m~) = (-) 



q s 7r 

2mr 




Ami, 



M(x)(8t R ) 



2:1 



(120) 

Here x is defined as x = m^/m 2 -. M{x) and N(x) are 
invariant dimensionless form factors that we have conve- 
niently normalized so that M(x), N(x) ~ 1 when x — > 1, 
see Eas . 1 1 54| and 1 1 5 51 of the appendix. When calculating 
the supersymmetric contributions to the asymmetry one 
has to take into account the renormalization of the Wil- 
son coefficients from the s usy scale down to the bottom 
mass scale. Following Ref. [88| we have included the NLO 
corrections using the generalized factorisation approach 
assuming that mr ~ 1 TeV. We note that all the Slr 
contributions to the low energy effective Wilson coeffi- 
cients arise originally from the Wilson coefficient C g (m~) 
in Eq. 11201 Therefore all the contributions to each effec- 
tive Wilson coefficient coming from the flavor violating 
soft susy breaking trilinear couplings (i.e. Slr. couplings) 
can be added up since they are proportional to the same 
gluino-squark form factor M(x). The resulting contribu- 
tion can always be written in the form, 



^SUSYI 

<t> K s Ilr+rl 



2 UkM(x) 
18m;,m~ 



UU32J 



(121) 

In our notation the coefficient f^x absorbs the depen- 
dency on the method of calculation of the hadronic ma- 
trix elements. We have calculated J^k using the gen- 
eralized factorisation approach, following Ref. jig], fyx 



is parametrized in the form f$K — fg^4 



is a factor 



associated with the momentum carried by the gluon in 
the corresponding penguin diagram, f q = mj/y' (q 2 ) ( in 
the rest of the paper we will assume that f q = \/2). X^, 
arises from the hadronic matrix elements. It is given by 
X<f, = 2F®^ K (m 2 ^) f 4,171 <p(pK •£</>)■ The numerical value of 
the parameter is irrelevant for our purposes because 
it cancels with the same factor coming from the SM con- 
tribution. The coefficient J^k could be calculated using 
other more recent and precise approaches which are avail- 
able in the literature: as the perturbative QCD approach 
[89| or the QCD factorisation approach [9(j • Nevertheless 
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it has been pointed out that in that case one would ob- 
tain sligthly different values for the relevant coefficients 
[zEOniOSIIIl Therefore for our purpose, which is to ob- 
tain a good estimate of the constraint on the new physics 
contributions, the generalized factorisation approach is 
precise enough. 

The contributions to the supersymmetric amplitude 
coming from the <5ll an d £>rr couplings can also be writ- 
ten in a similar fashion, 



1 SUSY I 



s I LL+RR 



47m? 



(0*Ll)23 + (Si 



KR)23) 



(122) 

L(x) is a dimensionless polynomial conveniently normal- 
ized such that L(x) — > when x — * and L(x) « 1 
when x — > 1. We note that in this case, since there are 
<5ll and $RR contributions coming from different Wilson 
coefficients, the coefficients of the polynomial L(x) de- 
pend on the method used to evaluate the hadronic ma- 
trix elements and to a lesser extent on the scale of the 
supersymmetric spectra but they do not depend on the 
flavor mixing structure in the susy breaking sector. We 
have evaluated the coefficients of L(x) numerically using 
the generalised factorisation approach following Rcf. 
We obtain approximately, 



L(x) 



GF 



cq + ci(x — 1), c ~ 1, ci « — v3, (123) 



in the limit x — > 1. Additionally, if one is interested in 
the limit x ~ 0, i.e. rrig <C rriq, L(x) is approximated by, 

L{x)gf = (— )4x (d + d\ lax) , rf ~ 18, d\ « 7, 

(124) 

where the coefficients co,i and do,i shown here are good 
approximations to the actual values calculated numeri- 
cally. The expressions in Eas. I12ll and 11221 are practical 
expressions of general interest irrespective of the form of 
the matrices Slr, Sll and <5rr. We see from Eqs. 11211 
and 11221 that one naively expects that the <5lr contri- 
butions dominate since the Sll and 5rr contributions 
receive in general an additional suppression factor of the 
order mf,/(5m~). Nevertheless for the model under con- 
sideration we obtain, 



(O23 + (6i R )* 32 = 4A(c^-^)^ C/3 , (125) 

b 



(SD23 + (4r) 23 = -4e»p'\. (126) 

We note that the total Slr contribution to S^Ks ^ s zero 
since XmKS^)^ + (^lr)32] = 0. Therefore in our model 
we find that only (LL+RR) couplings contribute to Sqk s ■ 
We find the following simple expression for the ratio of 
the dominant supersymmetric contribution to the ampli- 
tude over the Standard Model contribution, 



R, 



^Ill+rr 




L(x)\(5t 



LlJ23 



RRJ23\ 



(127) 



where 77^ is a coefficient independent of the supersym- 
metric parameter space given by, 



V2a 2 s 



^G F \V; b V ts \h^ 



(189 GeVf 



(128) 



Here h$ parametrizes the dependence of the SM contri- 
bution on the Wilson coefficients and hadronic matrix 
elements. We used the value for calculated numeri- 
cally in Ref. using the generalised factorization ap- 
proach (GF) . For instance if rriq — 500 GeV we obtain 
I^Ill+rr ~ 0.14L(x) \{5 d LL ) 2 3 + (Sinhsl which agrees 
with previous numerical calculations [741 . Isfjj . We note 
that Eq. 11271 provides some analytical insight in the de- 
pendency of the supersymmetric contributions on the su- 
persymmetric spectra, especially on the gluino squark 
mass ratio through the form factor L(x). 

Finally we will use the expression for ((Sf L )23 + 
(5^)23) in our model given in Eq. 11261 and the gen- 
eral expression for the amplitude A^Ks given in Eq. 11271 
to rewrite the constraint from Eq. 11191 Using for A, 
A k, 0.22, we obtain. 



212 GeV 



L(x) < 0.40 ± 0.26 (129) 



We note that the phase </>np as well as the strong phases 
difference S are not constrained by the data on quark 
masses and mixings. If 0np = this contribution to the 
asymmetry Sqk would cancel. Let us assume in the worst 
case scenario that </>np = 7r/2, 6 = ir/2 and p' = p = 9 
(which is the value for the large tan (3 scenario analyzed 
in Sec. IIII AT) . We would obtain only a mild constraint 
on the squark mass scale, of the order m~ > 1 TeV. 

Finally we would like to mention that, as it has been 
pointed out, in the case when the <5lr contribution is 
much smaller than the <5ll or 5rr contributions, the 
chargino contributions to the amplitude may be relevant 
since they could be of the same order than the 5ll and 
Srr gluino contributions [73, [ZS [7!J Hfl |8M H|| • A more 
precise calculation would require the inclusion of these 
contributions. 



C. CP asymmetry in B —* rj'Ks 

Recent results on the measurements of the CP asym- 
metries on the b — > s processes have reported possi- 
ble anomalies not only in B — » 4>K S but also in other 
processes, including B — ► T]'K S . The latest results 
from BELLE [Hi and BABAR [f| collaborations for 
the time dependent CP asymmetry coefficient S v ik 3 , 

(^i LLE ) 2 = 4 +0.06 ± 0.42 and (S^ AK ) ^ +0.50 ± 
0.32, seems to differ from the SM expectation. Combin- 
ing the results from both experiments one obtains the 
world average, (^abar+bellE) 2 0= 04 +Q 41 ± n p 

This has motivated the recent interest in the supersym- 
metric contributions to the CP asymmetry in the decay 
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B — ► r]'K s versus B — > </>if s [8a, |9j| as well as in cor- 
relations with other supersymmetric processes [85l l92j . 
It is known that because vector mesons (4>,p,- ■ •) and 
pseudoescalar mesons (71", K,rf , • • •) have opposite parity 
the B decays to these two final states will be sensitive to 
different combinations of the relevant Wilson coefficients 
93] . For instance, in supersymmetric theories the gluino 
loop effects coming from i5lr couplings will contribute by 
a factor proportional to ((<5lr)23 + (^lr)^) in the vec- 
tor case and to ((<5lr)23 — {Sim)^) in the pseudoscalar 
case respectively. For the model under consideration, 
the contributions from <5lr couplings, which exactly can- 
cel for the S<pK a asymmetry, not only do not cancel but 
dominate the CP asymmetry in the decay B —> rj'K s . 
We can obtain an expression for the $lr contribution 
to A^'Ks similar to A^k s in Eq. 11211 with the change 

((<$lr)23 + (<5lr)3 2 ) ((<5lr)23 - (Slr)* 32 )- Using the re- 
sulting expression we obtain the following simple formula 
for the supersymmetric contribution to S^k, , 




Using the values of A and lo determined from quark 
masses and mixings, and assuming that A b ps m~ we 
obtain the following constraint on S^'k, 



low-energy effective weak Hamiltonian [95 



1 



\C?\ 



+a g7 Im[C£Cf' 



[a 27 Im[C2(C 7 L * + C 7 fl *)] 



+ a 2g Im[C 2 (Cf+Cf )]] 



(132) 



where Cf = Cf ctt (m b ), = C^ ctt (m b ) and C 2 

multiply the chromo-magnetic dipole operators, O7 = 



167T 



Oa 



9s 

KiTT-' 



rSL&tivG^bft, and the 
current-current operator, 2 = SLJuqLQLj^bjj, respec- 
tively, and are the corresponding coefficients of 
the non-standard dipole operators involving right-handed 
light-quark fields, which appear in supersymmetric theo- 
ries. We will use the numerical values of the coefficients 
ay as computed using the parton model in Ref. [93| : 
a 27 w 0.0123, a g7 « -0.0952 and a 2g ps 0.0010." In 
order to explore the implications of supersymmetric fla- 
vor models it is useful to express the effective coefficients 
in terms of the new physics contributions to the Wil- 
son coefficients at the scale mw- To this end numerical 
expressions were given in Ref. |93|| including NLO renor- 
malization effects from mw down to the mj, mass scale, 

C 7 = C? + 7777 C 7 (m w ) + mg C g (m w ), (133) 
C g = C?+tfc C g (m w ). (134) 



\S<j > CpC 2< j 11 Ss\ 



1 TeV 



M(x) < 0.33 ±0.16. (131) 



We note that the phase </> as well as i5, the difference 
between strong phases, are not constrained by the data. 
If <fi pa this contribution to the asymmetry would cancel. 
In the worst case scenario, assuming that x ~ 1, <f> = tt/2 
and S = tt/2 the constraint depends strongly on the value 
of tan/3. For large tan/3, tan/3 = 50, we would obtain 
a mild lower constraint on the squark mass scale, m~^ > 
250 GeV. 



D. CP asymmetry in B — » X a ^y 



CLEO collaboration has set a range on the direct CP 
asymmetry in the b — > 57 decay, A^ 61 , at 90 % C.L. 
as A^ S7 = (-3.5 ± 13.5)% 25] while the BELLE col- 
laboration also set a range as 

^CP S7 = (-°- 8 ± 10 - 7 ) % 
[9^|. According to the SM theoretical prediction A^^ 

is smaller than 1% j^]. Therefore A^ 87 is an observable 
potentially sensitive to the presence of new physics. Fur- 
thermore it is expected that the experimental uncertain- 
ity will be reduced to less than 1% at a super B factory. 
^Cp S7 i n supersymmetric theories has received consider- 



able interest recently |2a, 123,122) 98] . It is known that a 
CP violating phase in the entries (<^lr,rl)23 or (Sf T] 00)23 
will generate CP violation in the decay B — > X s j |9al99| . 
The direct CP asymmetry in b — > sj decay can be writ- 
ten in terms of the effective Wilson coefficients of the 



Here 7777 = 0.67, r\ 7g = 0.09 and i] g = 0.70. The super- 
symmetric contributions to C g were given in Eq. 11201 



Cf(m 7 ) 



2m~? 

b 



b 



M 3 (x)(5t L ) 



23 



L Ml (x)(Sl R ) 23 



(135) 

where Mi(x) and M 3 (x) are dimensionless form factors 
defined in Eqs. 11581 and 11591 of the appendix. M\{x) 
and M 3 {x) have been normalized to 1 when x — > 1. 
We note that for simplicity we have defined C 7 as the 
whole coefficient accompanying the operator 7 . There- 
fore in our notation the SM contribution to the Wil- 
son coefficient C 7 at mw is given by C 7 M (mw) — 
~^/2m b GFV t * s VtbK (xt), where xt — mt/mw and K(x) 
is a dimensionless form factor given in Eq. 11601 of the 
appendix. The supersymmetric contributions to C 2 are 
neglibible. We will use the SM value, C 2 {m b ) ps 1.11 x 
GFV t * s Vtb/V%- It is straigthforward to obtain and Cf 
by the exchange L «-> R in the expressions for C 7 and 
Cg. Barring cancellations between <5lr and 5ll terms 
we will obtain an approximate bound from the LR con- 
tribution. We can see from Eq. 11321 that the total <5lr 
contribution is proportional to a coupling of the form 
((<5lr)23 + (^LiOJb)- We obtain the following approxi- 
mate expression for the asymmetry, 



.4 



b — >S7 
CP 



LR+RL 



(-) 



a s 7r C 2 lm[(S hR )* 23 + (5 hR )* 32 JA(x) 



2m~ 

b 



|cf M | 



(136) 
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For the 5ll,r.r couplings we will obtain a similar ex- 
pression proportional to the coupling ((<5ll)23 + (^RjOsb)- 
Here A(x) is a dimensionless form factor defined by, 



mixing couplings contribute to the relevant Wilson coef- 
ficients. The b — > S7 decay rate is also proportional to 
the CV Wilson coefficients, 



A(x) = {a 27 rm-M 1 {x) + 

a2g(V7g^M 1 ( X ) + T} g M(x))) 



For the model under consideration, 



vA h 



(O23 + O&OSa = (-)4A s ^^ C/3 . 



(137) 



(138) 



The SM contribution to the effective coefficient, Cj is 
related with the Wilson coefficient at the mw scale by 
a renormalization factor, — ?7&wCf M , which can be 
extracted from Ref. 93]. Assuming that x ~ 1, i.e. m— ~ 
mr, and using the values of A and u> as determined from 
quark masses and mixings we obtain the constraint, 



LR+RL 



S<j>Cp 




<0.1 (139) 



In the worst case scenario, assuming that w 1, 4j « 
mr, tan/3 1 the current experimental bound requires 
to- > 230 GeV. On the other hand, for large tan/3 one 
would obtain a milder constraint. We would like to point 
out that the phase <f> is not constrained by the CKM 
phase. If Sd, <C 1 the squark masses would not be con- 



b — >S7 
CP ■ 



One would naively expect that the 

b — >S7 



strainted by A 

<5lr gives the dominant contribution to A° c ^, s ' 1 because 
of the toj, jv suppression factor of the (5ll contributions 
to C7 and C g . Nevertheless for the model under consid- 
eration the contributions coming from <5ll.rr couplings 
are of the same order of magnitude. For the model under 
consideration, 



(SD23 + (O32 = -4Ap'e" 
We obtain a similar expression, 



(140) 



r(6- 



ST ) cx (|Cf| 



+ \c?\y 



(142) 



The current world average of the CLEO |102j and BELLE 
[l03] results is given by B(b -> S7) cxp = (3.3±0.4) x 10~ 4 , 
which can perfectly be accounted for the SM theor etica l 
prediction, Bib -► S7) th co = (3.29 ± 0.33) x 10" 4 [l0^ . 
which leaves a small window open for new physics. There 
is no SM contribution to Cy . A full expression for the 
main supersymmctric contributions, i.e. gluino exchange, 
to this branching ratio were first given in Ref. jlOO j . Fur- 
ther improvements in the calculation as chargino dia- 
grams and QCD corrections were subsequently included 
|101| . Therefore if the supersymmetric contribution is 
just a correction to the SM one we can expand in powers 
Cj susy /Cy M | and obtain to leading order, 



of R SJ = 



T{b -* s-y) 
T(b -> s 7 ) SM 



- 1 w 2R 



«7- 



(143) 



Allowing the supersymmetric contribution to saturate 
the 2a experimental uncertainity we obtain to leading 



order in R SJ the constraint, 



R < 



A(Br(6^ S7 )l eX p) 
Br(6^ 57) I 



£ 0.12, 



(144) 



cxp 



Using the expression for the supersymmetric contribution 
to Cf from Eq. 11351 we obtain for the i5lr contribution 
to R sl the expression, 



R 



S7ILR 



(S d LR ) 23 \M 3 (x). 



(145) 



Here rj sl is a coefficient independent of the supersym- 
metric parameter space with mass units. Using the SM 
expression for Cf and the expression for (<*>lr)23 that our 
model predicts for x ~ 1 , 



^CP 



<0.1. (14!) 

LL+RR \ TOr I 



The constraint on the squark spectra depends on the 
value of p' . For the large tan /3 case we noted in Sec. lIII Al 
that p = to 2 /to- 9. If p' — p we would obtain the con- 
straint to- > 300 GeV. This constraint could be avoided 

b 

if 4> ~ or rf < 1. 



e. r(6 -> si) 

The supersymmetric contributions to the b — > S7 de- 
cay are indirectly correlated with the CP asymmetries 
in B — > 4>K S and b — > S7 decays since the same flavor 



\ I ^ OA U m 

|(0lr)23| « 2AC/3 — — 
we obtain the constraint, 



TOr to- 




c /3 



130 GeV 



< 0.12. 



(146) 



For the large tan 3 case with to— ss mr examined in 

_ 9 b 

Sec. lIII Al . tan 3 = 50, m was required to be m ~ 3mr. In 
this case we obtain the following constraint on the squark 
mass scale, m~ > 370 GeV. An analysis with similar re- 
sults can be implemented for the <5(ll+rr) contribution 
to T(b — > sj). In this case, 



23 



2Xp'. 
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We obtain the constraint, 
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(147) 



mr examined in 

b 



For the large tan/? case with m~ 
Sec. MI Al , tan/3 = 50, p was required to be p » 3. 
\i p' = p we would obtain the following lower bound on 
the squark mass scale, > 140 GeV. 



VIII. CONCLUSIONS 

We have shown that generic supersymmetric flavor 
models exist for the radiative generation of fcrmion 
masses, mixings and CP phases. We have studied in 
detail the phenomenological implications of a particular 
supersymmetric flavor model for the radiative generation 
of first and second generation quark masses, focusing our 
attention especially in the CP violating phenomenology 
The basic idea underlying this kind of flavor models is 
that the flavor breaking fields are also supersymmetry 
breaking fields. 

We have shown that these models generically solve the 
SUSY CP problem in a very simple fashion. The one-loop 
contributions to EDMs are exactly zero due to the real 
character of the relevant parameters while the two-loop 
contributions are suppressed. 

Our main goal was to present a flavor model as predic- 
tive as possible. To this end we have proposed a partic- 
ular implementation of this scenario using a U(2) flavor 
symmetry where the required hierarchy of flavor breaking 
vevs is expressed in powers of a single parameter, A. As 
a consequence the model generates quark Yukawa matri- 
ces that contain only three parameters, A, 9, 7 and can fit 
the data with precision. Therefore the quark masses and 
mixings determine the amount of flavor violation in the 
soft sector requiring a very heavy susy espectra especially 
to overcome the constraints on Am^c , e and Am^. 

We would like to emphasize that this case study can be 
considered the worst case scenario from the point of view 
of FCNC constraints. Between the extreme case study 
here considered and the usual models with scalar flavor 
breaking vevs there is a continua of possibilities which 
would ameliorate the FCNC constraints. For instance, 
we could increase the number of parameters in the fla- 
vor breaking sector, use a different flavor symmetry or 
generate radiatively only the first generation of fermion 
masses. If that was the case one could lower considerably 
the constraints on the sfermion spectra, probably at the 
price of decreasing the predictivity of the flavor model. 

We believe these models are an scenario worth of a 
more dedicated exploration. They generically allow us 
to simplify the "flavor vacuum" , or in other words the 
hierarchies of the flavor breaking vevs, through the intro- 
duction of a natural hierarchy, the loop factor, and they 
offer a new insight in the SUSY CP and flavor problems. 



For completeness we include expressions for the dimen- 
sionlcss form factors that were used in the main text. The 
form factor F(x, y, z) is defined as, 



F(x,y,z) 



(x 2 y 2 In + V 2 z 2 In p + z 2 x 2 In §3 



(x 2 — y 2 )(y 2 — z 2 )(z 2 — x 2 ) 



f(x), g[x) and h(x) are given by, 

'(x 2 - 82;- 17) 6(l + 3z)ln(a;)) 



(148) 



m = — 



(X -l)* 



(x-iy 



g(x) = 10 
h(x) 



(.t 2 + 10.t + 1) 6x(l + x)\n(x)) 



(X -l)* 
llg(x) - 6/(a;) 



(x - 1) E 



(,149) 
, (150) 
(151) 



These functions appear in the calculation of the super- 
symmetric contributions to the Wilson coefficients The 
functions f(x), g{x) and h{x) have been conveniently 
normalized so that in the limit x — > 1 they tend to 1. 
Approximate expressions in the limits x — > 0, 1 are given 
by, 



/(*) = 



_ J i-±c + 0(C 2 ), x^l,{C = x-l) 



lte "(17 + 61n(a;)) +0(x 2 ), x -> 0, (152) 



1-C + 0(C 2 ), a: -» 1, (C = a: - 1) 
10(1 + 2x(3 ln(ir) + 7) + 0(x 2 ), x -► 0, 



(153) 

N(x) and M(x) are dimensionless form factors given by, 
~{x 2 + I72x+ 19) 



N(x) 



M(x) 



36(i- l) 4 

(x 2 - 15x - 18) 



-x ln(x 



6(x - l) 5 
(54x 4 + 216a; 3 - 287x 2 - 8x + 1 

9 (a: - l) 4 
2,„^( 36a;2 ~ 19a; -21) 



(154) 



-2*M«)- 3(x _ 1)5 
In the limit x ~ 1 M(x) is given by, 



(155) 



31 233 

M{x) = a Q + ai {x-\)+0((x-l) 2 ), a = -,«i = — . 

(156) 

If one is interested in the limit x ~ 0, i.e. nig <§C rrig, it 
is also possible to obtain an approximate expression for 
M(x), 

M(x) = -Vi (6 + hx + 0(x 2 ), ) b = 1, 61 = -4. 
9 

(157) 
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The functions M\{x) and M 3 (x) are defined by, 
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M 1 {x) 



M 3 (x) 



12x 2 ln(x) 6(1 + 5x) 



{1-xY 



+ 



10^/x 



(1-8 



(1 - x) 3 
x - I7x 2 ) 



(158) 



ln(x; 



(18x 2 + 6x 3 ) 



(x - 1) E 



159) 



Here Mi{x) and M 3 (x) have also been "normalized" so 
that in the limit x — * 1 they tend to 1. Finally the 
dimensionless form factor K(x), which appears in the SM 
contribution to the Wilson coefficient CV, is given by, 



K{x) 



2(x - 1) 



(8a; 2 + 5x - 7) 
12 



2X 2 
3 



ln(x) 

(x~^T) 
(160) 
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